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g^' Abstract 



It is shown that in the Einstein- Yang-Mills (EYM) theory, as well as in the 
pure flat space Yang-Mills (YM) theory, there always exists an opportunity 
5^ i to pass over the potential barrier separating homotopically distinct vacuum 

sectors, because the barrier height may be arbitrarily small. However, at 
low energies all the overbarrier histories are suppressed by the destructive 
interference. In the pure YM theory the situation remains the same for 
any energies. In the EYM theory on the other hand, when the energy is 
large and exceeds the ground state EYM sphaleron mass, the constructive 
interference occurs instead. This means that in the extreme high energy 
limit the exponential suppression of the fermion number violation in pure 
YM theory is removed due to gravitational effects. 



^On leave from Physical- Technical Institute of the Academy of Sciences of Russia, 
Kazan 420029, Russia 



1 Introduction 

The discovery of the vacuum periodicity in a Yang-Mills (YM) theory [|l| 
gave rise to the theoretical description of baryon number non-conservation 
[0] caused by the anomalous violation of chirality 0. The standard approach 
to understanding this phenomenon uses semiclassical barrier penetration be- 
tween topologically inequivalent vacua, where the tunneling solutions are 
Euclidean instantons [Q. The amplitude of such a penetration at zero en- 
ergy is of the order of exp{—S^'^'^''), which is incredibly small because of the 
large value of the instanton action S^^'^K In a pure Yang- Mills theory one 
can not expect the amplitude to become large with growing energy, because 
the theory does not contain any dimensional parameters to set the energy 
scale. 

When scalar fields are involved, as in Yang-Mills-Higgs (YMH) theory, a 
natural energy scale arises, and one may think that the transition amplitude 
will not be small if the energy is of the order of some threshold value. In- 
deed, the existence of the sphaleron solution in the electroweak YMH theory 
[^] implies that the fermion number violating processes may become unsup- 
pressed if the energy (or temperature) is of the order of the sphaleron mass 
[^, [7|. The idea is that, instead of tunnel penetration under the barrier, the 
field configuration can pass over it and the value of the sphaleron energy is 
the minimal barrier height. 

If the energy is very high, then the symmetry is restored. The YMH 
sphalerons then dissociate, and the fermion number violating processes again 
become suppressed |]^. However, another universal energy scale arises nat- 
urally at extremely high energies, namely, the gravitational Planck scale. 
The existence of this new scale leads one to expect that the fermion number 
violation in the YM theory will become significant when the energy is com- 
parable with Planck's energy. This expectation is supported by the discovery 
of particle-like solutions in the Einstein- Yang-Mills (EYM) theory 0, which 
have been interpreted as sphalerons |T^. One might think that the typical 
mass of these EYM sphalerons, similarly to the situation in the YMH theory, 
determines the transition threshold. 

It turns out however, that, contrary to the situation in the YMH theory, 
none of the EYM sphalerons relate to the highest energy point on a minimal- 
energy path connecting inequivalent EYM vacua, for the potential barrier 
between distinct vacua in the EYM theory can be arbitrarily low [T^. Indeed, 

1 



let us consider first the fiat space YM theory and define a one-parameter 
continuous family of static field configurations, A{x, A), such that when the 
parameter, A, goes from zero to unity, the field A{x, A) interpolates between 
two vacua with distinct winding numbers. We will call such a family of 
fields a vacuum-to-vacuum (VTV) path. The simplest choice is ^ A{x, A) = 
XA^^\x), where A^^\x) is a pure gauge with unit winding number. The 
energy for this family, U{X), has the typical barrier behaviour as it vanishes 
for the vacuum values, A = 0, 1, and reaches a maximum in the middle at A = 
1/2. This allows one to say that distinct vacua are divided by a barrier. Next 
note, that continuous scahng deformation of the family, A{x, A) -^ A{Px, A), 
alters the energy as follows: U{X) — >■ f3U{X), and hence the barrier height, 
(3U{l/2), can be made arbitrarily low by a proper choice of the scaling factor 
p. It is clear that the inclusion of gravity can only reduce the energy and 
hence does not change the situation (a consistent way to determine of the 
gravitating energy demands solving the initial value constraints; see below). 

Thus, the EYM sphalerons, having finite mass, are not related to the 
mmma/ barrier height. One may wonder then what their physical meaning is. 
Indeed, as the barrier height may be arbitrarily small, an overbarrier passage 
is possible at any low energy and not only at the (ultra-high) sphaleron 
energy. In other words, there is no need to climb the mountain in order to 
reach the neighbouring valley as it is possible to go around it. This certainly 
obscures the threshold interpretation of the sphalerons. 

However, there exists also a counter example to this objection. Indeed, 
in the fiat space limit, where the overbarrier passage is also possible at any 
non-zero energy, the transition amplitude nevertheless always remains small 
(in the pure YM theory), as all fermion- violating Green functions contain 
the small instanton factor, regardless of the value of the energy 0. There- 
fore, one may infer that the possibility to pass over the barrier does not yet 
guarantee in itself an enhancement of the transition rate. To understand why 
this may happen one has to consider not only one path interpolating between 
distinct vacua, but all such paths, which corresponds to the inclusion of all 
the possible transition histories. The contribution of a single path to the 
transition amplitude may not be small provided that this path lies entirely 
in the classically allowed region, that is, the energy exceeds the barrier height 
for the given path - we shall call such a path an overbarrier path. However, 
when one takes the sum over all histories, it may happen that the contri- 
butions of different paths cancel each other. The latter occurs when the set 



of all paths does not include a path whose contribution into the transition 
amplitude is stationary with respect to small variations of the path itself. 
Such a situation is encountered in the flat space pure YM theory, where all 
overbarrier histories are strongly suppressed by the destructive interference, 
regardless of the value of the energy, and a stationary phase path, if any 
exists, may only be underbarrier, so that the total transition amplitude is 
always small. 

One may infer from this example that for the EYM theory all the fermion 
number violating transitions should also be suppressed in the low energy 
limit, despite the existence of the overbarrier paths, for gravitation can not 
essentially alter the dynamics at low energies, and the theory should resemble 
the fiat space YM theory in this limit. Therefore, the vanishing of the great- 
est lower bound of the barrier height may be not so fatal for the threshold 
interpretation of the EYM sphalerons. Indeed, as we will see below, the re- 
semblance between the YM and the EYM theories only holds at low energies. 
In the high energy limit, the situation changes. Firstly, at high energies the 
set of the overbarrier paths in the EYM case becomes larger in comparison 
with that appearing in the fiat space theory, because the gravitational bind- 
ing tends to reduce the barrier height. In addition, and this is the crucial 
point, when the energy exceeds the ground state sphaleron mass, gravity in- 
troduces a stationary phase overbarrier history. This leads us to conclusion 
that the transition becomes unsuppressed at high energies and restores the 
initial threshold interpretation of the sphalerons. 

The main goal of this paper is to argue, using the quasi-classical argu- 
ments, that the suppression of the fermion number violating transitions in 
the EYM theory is removed when the energy is larger than the ground state 
sphaleron mass (see also |3^). The theoretical tool relevant for our purposes 
has been developed by Bitar and Chang in their real time analysis of tunnel- 



ing transitions in the fiat space YM theory |11]. The essence of this method 
is to introduce firstly a VTV path in the configuration function space of a 
theory. Then one assumes that the transition between distinct vacuum sec- 
tors of the theory proceeds along this single VTV path. This reduces the 
problem to a one-dimensional quantum-mechanical task and allows one to 
write down immediately the corresponding WKB amplitude. The last step 
is to take the sum over all VTV paths. 

We apply Bitar and Chang's technique to represent (in the lowest order 
WKB approximation) the amplitude of the fermion number violating tran- 



sition in the EYM theory at non-zero energy (a similar idea was exploited 
also within the context of the YMH theory; see Ref. |T^). Following this line, 



we introduce a set of VTV paths connecting homotopically distinct vacua in 
the EYM configuration function space (see Eqs.p3|),(^) below). Our basic 
approximation is the assumption that the quantum transition between dis- 
tinct vacuum sectors occurs only along these paths. We estimate the partial 
WKB amplitudes for the paths introduced and specify those paths whose 
amplitudes are not small. Such paths lie entirely in the classically allowed 
region and they relate to overbarrier passages. Remaining paths pass through 
the classically forbidden region and we exclude them from our consideration 
because the corresponding partial amplitudes are small due to the tunnel- 
ing suppression. The total transition amplitude is the sum over all partial 
amplitudes which can be estimated by making use of the stationary phase 
approximation. We show that the overbarrier paths always exist, that is, 
the overbarrier passages are possible at any non-zero energy. However, when 
the energy is small, the set of the overbarrier paths does not contain a sta- 
tionary phase path. This means, that the sum over all overbarrier histories, 
being evaluated via the stationary phase approximation, will be small, as the 
stationary point is absent, so the overbarrier passages are suppressed by the 
destructive interference. On the other hand, when the energy is large and 
exceeds the ground state sphaleron mass, a stationary phase path appears 
on the set of the overbarrier paths. This allows us to conclude that, when 
the energy is large, the total transition amplitude in the EYM theory is not 
small (because the stationary point is present). Thus, the suppression of the 
fermion number violation in the EYM theory is removed when the energy is 
of the order of the EYM sphaleron mass. 

The rest of the paper is organized as follows. 

In Sec. II we investigate the homotopical classification of the vacua arising 
in the SU(2) EYM theory and estimate instanton contributions to the fermion 
number violation rate. In Sec. Ill we introduce a set of vacuum-to-vacuum 
paths in the EYM configuration function space such that the time evolution 
along these paths leads to a change of the topological winding number. In 
Sec. IV we apply Bitar and Chang's technique to obtain the WKB transition 
amplitude. To select the overbarrier paths we study in Sec.V the structure 
of the potential barrier dividing vacua in the EYM theory. In the vicinity of 
the n-th sphaleron solution, we find that locally, the potential barrier surface 
is a saddle with one transversal and n longitudinal negative directions. Ex- 



ploiting local properties of the potential barrier surface we show in Sec. VI the 
existence of the (approximative) overbarrier stationary phase path when the 
energy exceeds the ground state sphaleron mass. Some concluding remarks 
are made in Sec. VII. There are also three appendices included. Appendix A 
contains an explanation of our choice of the gravitational action, and also 
details of the computation of the action for the spherically symmetric case. 
The Appendices B and C include various formulae relevant for the descrip- 
tion of the spherically symmetric EYM fields and an explicit procedure for 
the first and second variation of the ADM energy functional. 

Our sign conventions used are those of Landau & Lifshitz ||37|| . 



2 The EYM vacua and the instanton esti- 
mates 

The starting point of our considerations is the notion of the topological vacua 
of the EYM fields. The topological vacua of the YM field have been first 
introduced in Ref. [|l| within the context of the fiat space theory. It turns 
out that inclusion of gravity modifies the analysis in |jl| only slightly, as long 
as one does not take into account the topological effects of the gravitational 
field itself. 

Consider the action of the EYM theory with the SU(2) gauge group 

Seym = Sq + Sym, (1) 

with the gravitational part 

^<= = - lib / ""^"'^ - lib £(*'"r"^ - ^'-'r-^'''^- P' 

and the Yang-Mills contribution 

Sym = -7^ f trF,,F^^^^gd''x. (3) 



2^^ 

Here G is Newton's constant, g is the gauge coupling constant, F^^ = d^A^ — 
dyA^ — i[A^, Ay\ is the matrix valued gauge field tensor, A^ = y4^r"/2, and 
r" (a = 1,2,3) are the Pauli matrices. We use the gravitational action. 



which does not contain second derivatives of the metric; for explanations of 
our choice see Appendix A. 

The action (|I|) has many stationary points, including solutions with non- 
trivial space-time topology such as black holes ||15|. Consider a sector with 



spacetime manifolds carrying the trivial R^ topology. Define vacua in this 
sector as the stationary points with zero ADM energy. As follows from the 
positive energy theorem [|l^, such vacua must have a flat metric and hence 
the YM field is a pure gauge. In the Aq = gauge one can represent the 
vacuum fields as 

Aj{x) = lUdjU'^, g^,y{x) = ri^y, (4) 

where rj^y is the flat space metric, U{x) is a SU{2) valued function, and 
spatial index j runs over 1,2,3 [x = x^). 

Such vacua are to be the boundary conditions for the path integral. It 
seems that not all of them are physically important, but only those satisfy- 
ing some additional topological conditions, which have the origin in instanton 
physics. Following the standard line [|ll, we pass to the Euclidean sector and 
use the instanton arguments to clarify the structure of the EYM vacua. One 
should note that restricting ourselves only to spacetime manifolds with the 
R^ topology, we exclude from consideration all non-trivial (i.e. carrying non- 
zero topological charges) gravitational instantons and vacua. A complete 
investigation of quantum gravity effects requires inclusion of instantons and 
vacua with non-trivial topological indices, both for the gauge and gravita- 
tional fields. This lies, however, beyond the scope of our present analysis 
(description of the non-trivial gravitational vacua has been done in Ref.||17||). 

Consider an instanton interpolating between two EYM vacua. The finite- 
ness of the instanton action implies the finiteness of both the matter and the 
gravitational parts of the action as both of them are positive. Indeed, the 
matter action is positive. The gravitational Euclidean action may be arbi- 
trary, however the positive action conjecture implies that it should be positive 
on shell, where the condition i? = holds |TB[. Finiteness of the matter part 



of the action imphes that the YM field tends to a pure gauge, iUdjU , at 
Euclidean infinity, S'^^^i, which defines a map 5*1;^^; — > Sgm^2)y ^s the SU{2) 
group is also a three sphere 0] . All maps S^ —>■ S^ fall into homotopy classes 
labeled by an integer degree of the map. In the case under consideration the 



degree of the map is the Pontryagin index 

y=^^jtrF,,F^^^Sx, (5) 

where the dual tensor is F^^ = ^1=6^''°'^ Fap {e^^"^^ = !)• It is worth noting 
that this equation in fact does not involve the metric, so we infer that instan- 
tons interpolating between two vacua in EYM theory may be still described 
in terms of gauge invariant, topological indices of the YM field. 
Let us introduce the Chern-Simons current 

Stt^ v^ 3 

whose divergence is 

V^i^" = ^trF^,F'^^ (6) 

where Vq, is the covariant derivative with respect to the spacetime metric. 
Pass to the Aq = gauge, which violates the equivalence between the Eu- 
chdean coordinates. Using Eqs.(|5|),(||), one obtains 



^K^d^x + dtj KdT. = / ^K^d^x 



=+oo 



= k{oo) — k{—oo). 

-oo 

(7) 
The surface integral in this expression can be omitted as the YM field tends 
to a pure gauge at the spatial infinity. Indeed, for a pure gauge field taken in 
the Aq = gauge, K^ is the only non-zero component of the Chern-Simons 
current. So, the surface term in Eq.(0) vanishes when the boundary moves 
to infinity. In the limit r -^ ±cx3 the gauge field tends to a pure gauge, 
iUdjU^^, which allows one to write 

k= f K^d^x = -^ j e'^hrUdiU-^UdjU-^UdkU-^d^x. (8) 

Thus, one can attach to any EYM vacuum a number, k - the Chern-Simons 
charge - and the difference of two such numbers for vacua connected via an 
instanton is always an integer and gauge invariant as well. The value of k 
itself is neither gauge invariant nor an integer, however if the function U{x) 
in Eq.(^ satisfies the following additional condition [|l| 

lim U{x) = 1, (9) 

|a?|— >oo 



then k will be an integer and may be treated as the topological winding 
number. Indeed, the condition (^) implies that U{x) may be viewed as a map 
S spatial -^ S^su(2)^ ^hcrc ^,%„« is the compactified 3-dimensional space ^ 
|21| , so an integer degree of the map again appears, which is exactly given by 
Eq.(|^). Also, k will be invariant with respect to small gauge transformations 
generated by functions obeying (^ and having zero Chern-Simons index (P). 
Notice, that the condition (^ arises naturally. Indeed, for the trivial 
vacuum one has = Aj = iUdjU~^, which implies that U{x) = 1 (up to a 
global gauge transformation), so the condition @ holds. During the time 
evolution one has in the Aq = gauge 

= Ao{t,x) = lim Ao{t,x) = \im iUdtU-\ (10) 

|x|— >oo |x|^oo 



which gives |^ 

lim U{t,x) = 1. (11) 

Thus, the condition @ holds for any vacuum which may be connected with 
the trivial vacuum via an instanton. 

So, one can see that the vacuum classification in the EYM theory re- 
sembles that arising in the flat space YM theory. The EYM vacua can be 
classified in terms of the integer winding numbers of the gauge field. Grav- 
itation does not introduce any major peculiarities as long as one considers 
only manifolds with i?^-topology. 

It is natural to ask how large the amplitude of transition between two 
EYM vacua with distinct winding numbers is. As is known, in flat spacetime 
this amplitude is small because the corresponding instanton action is large 
[Q]. It is easy to see that the inclusion of gravity does not improve the 
situation. Indeed, the inequality 

tr{F^, ± F^,Y = 2tr{F^,F^'' ± F^^F^ > 

implies that the total EYM action satisfies the following relation 

SEym = A jtrF^^F^^^d'x + ^1"^' > ^kl + ^l"'^'. (12) 

As the Euclidean gravitational action is positive on shell, it is clear that 
the inclusion of gravity may only lead to an additional suppression of the 
instanton transitions between sectors with distinct winding numbers. 



However, as was first argued within the context of the standard model 
[^], instanton estimates for fermion number violating amplitudes may break 
down in the classical, many-quanta limit, where the sphaleron estimates are 
relevant. One may reformulate these arguments in the EYM case as follows. 
Let quarks, q and leptons, /, be involved. Then it follows from Eq.(|l2D that 
the fermion number violating amplitude, < qqql >, is always small. Indeed, 
such an amplitude may be obtained from the generating functional whose 
path integral always gives rise to the small factor exp{—S^^'^^). However, 
the inclusive amplitude, < qqqlA^h^ >, may not be small, where A and h 
stand for the gauge bosons and gravitons, and m and n are large numbers 
(see Ref.0 for details). Notice, that such an inclusive amplitude arises nat- 
urally in a process mediated by a sphaleron, as the latter, being a classical 
object, involves large number of quanta. In other words, one may expect 
that the collision of two high energy fermions produces the sphaleron as an 
intermediate state. Then it decays giving a large number of gauge and grav- 
itational quanta and also other fermions as a side effect due to the anomaly. 
A straightforward evaluation of the amplitude for such a process can be done 
with the use of perturbation techniques 0, but we shall proceed along a 
different line and apply instead non-perturbative approximation. 

3 An example of the real time winding num- 
ber changing history 

Now, we return to the Lorenzian sector and introduce a family of paths lying 
in the EYM configuration function space and connecting homotopically dis- 
tinct vacua. In the next section this will allow us to apply Bitar and Chang's 
technique in order to represent an amplitude of the winding number chang- 
ing transition. The considerations in this section follow the line sketched in 
Ref.|l|. 



We shall consider spherically symmetric fields and so we need the corre- 
sponding expressions for the field potentials. Let us use the following param- 
eterization of the spherically symmetric gravitational field 

ds^ = Rl{{l - —)a^de ^^ - r\d^^ + sinHd^^)}, (13) 

^ r 1 — 2m/ r 



where Rg = y/Anlpi/g is the natural length scale arising in the EYM theory 



with Ipi being Planck's length, all other quantities in Eq.(|l^) being dimen- 
sionless, and m and a being functions of the radial coordinate, r, and time, 
t, as well. Topological triviality of the space time manifold implies that the 
metric must be regular and asymptotically flat; this means that m, a are 
smooth functions satisfying the following boundary conditions 

o" ^ 1, m —>■ const as r ^ oo; m = o{r), a = ctq + o(r) as r — ;► 0. (14) 

The spherical YM field is given by Witten's ansatz p2|, which can be 
represented in the following form 

A = WqLi dt+WiLi dr+{p2L2~il-pi) L^} di3+{{l-pi) L2+P2 L-^jsin^ dip, 

(15) 
where 

Li = siwdcosip (— ) + sin^sirup (— ) + cosi) ("t)' 

L2 = d^Li, L3 = ——d^Li, (16) 

Wq,Wi,Pi,P2 are functions of r and t. It is convenient to use also another 
parameterization of this ansatz introducing functions Qq, Qi, f,a as follows: 

Wo = fio + d, Wi = Qi + a', pi = fcosa, p2 = fsina, (17) 

where dot and prime denote differentiation with respect to t and r corre- 
spondingly. One may see that only three of these function are essential, the 
fourth one, a, is actually a pure gauge parameter, we shall call it an angle 
parameter of the field. The gauge transformation 

A^ ^ U{A^ + td^)U-' (18) 

with 

U = exp{ip{t,r)Li} (19) 

preserves the form of the field (|15|) , (|T7|) leaving functions Qq, Qi, f invariant 



and altering the angle parameter as follows 

a -^ a + (3. (20) 

Also, gauge invariant components of the energy-momentum tensor of the field 
(|I5|) (see Eq.( |Bl[ )) do not include a. 
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If one puts in Eq.(|T5|),(|T^ fio = ^i = 0, / = 1, allowing for the a to be 
arbitrary, then the energy of the field vanishes, so the field will correspond 
to a pure gauge, iUd^U^^, where U is given by Eq.(|Typ with /3 replaced 
by a. The topological vacua arise if one imposes the following additional 
restrictions on the angle parameter: d = 0, which insures the Aq = gauge 
condition, and also 

a(0) = 0, a(oo) = -27rfc. (21) 



Indeed, the Chern-Simons charge for the field (|T5D (see Appendix B) is 

/\ poo 
^rz^K'^d^x = —j^ dr{W,ipl+pl-l)+p2p[ + {l-pi)p'2}; (22) 

in vacuum this expression takes the form 

-| oo 

— {sina — a) , 
2n 

which is the integer, k, provided that the condition ( PT| ) holds (the condition 
q;(0) = in Eq.(|2T|) insures regularity of the field at the origin). It is clear 
then that the function (|T^ with time-independent f3 vanishing both at the 
origin and infinity generates small gauge transformations as the condition 
( plD remains unchanged. 

Summarizing, we conclude that the gauge field (|T5|),(|1^ with fig = ^i = 
/ — 1 = and arbitrary time-independent a obeying the condition (|2T|) , 
represents the topological vacuum, iUdjU~^ with U = exp{iaLi), and the 
winding number is k. 

We consider now an explicit example of the winding number changing 
evolution of the gauge field [1^ : 

A^{t, ^) = ^(1 - K{r))Ud^U-\ U = exp(2^ALl), (23) 

where A is some function of time, A(to) = 0, A(ti) = vr, where to < ti, and 
K{r) satisfies 

K{0) = 1, A^(oo) = -1. (24) 

It is clear that the field (^) is spherical and hence it can be represented in 
the form (|l^) giving for the corresponding functions 

Wo = il-K)X, Wi = 0, pi = cos^X + Ksin'^X, p2 = {l-K)sinXcosX, (25) 
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and also 

_ ; K{K' - 1) _ K'ctgX 



/■2 • 2\r^2 , 2\ , {I — K)ctg\ 

f = sm XK + cos A, tga = — ^r-— . (26) 

•^ ^ K + ctg^X ^ ' 



One may see that the field (^31) vanishes in the beginning, A = 0, as well as at 
the end, A = tt, of the evolution (provided that A(to) = A(ti) = 0). Thus, at 
first sight, this field corresponds to a loop in the configuration function space 
interpolating between the same trivial vacuum. However, the crucial issue 
here is the right choice of the gauge. The homotopical vacuum classification 
is based on the Aq = gauge condition [Q , |2ll] , [^ , so one has to pass to this 



gauge. To do this, perform the gauge transformation pO| ) and demand 

Ao = (fio + a + /9)i^i = 0, (27) 

where VLq and a are given by Eq.(^). In this way we obtain the equation for 
the gauge parameter [3: 

f]o + a + /5 = A(l-ir)+/3 = 0, (28) 

which yields 

(5 = {K~1)X + B{t), (29) 

where the function B{r) may be put equal to zero using the residual freedom 
of the time-independent gauge transformations, which preserve the Aq = 
gauge condition. Thus, in the new gauge, the invariant functions VLq.VLi, f 
are still determined by Eg. (^61) , but the angle parameter a is replaced by 



ao = {K-l)X + arctg^\f^. (30) 

The corresponding expressions for the functions Wo,Wi,pi,p2 in the new 
gauge are 

Wo = 0, Wi = XK', pi = cosXcosKX + KsinXsinKX, 

P2 = cosXsinKX — KsinXcosKX. (31) 

One can see that the field ([15|) specified by these functions vanishes at A = 0. 
On the other hand, at the end of the evolution, at A = tt, the field (p!5|),(pl]) 
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is the non-trivial pure gauge. The angle parameter of the new vacuum can 
be defined from Eq. (|30D : 

ao{r) = 7c{K{r)-l). (32) 

One may see also from Eqs.(^),(^) that this indeed corresponds to a vacuum 
with unit winding number. 

Thus the field (El) indeed interpolates between distinct vacua. To show a 



change of the winding number during the evolution (p3D one has to estimate 
the Chern-Simons number as a function of A (see also [^] for discussion of 
the subject in the case when arbitrary gauge groups are involved). Inserting 
functions (|3T1) into Eg . (p2[) and using Eq.(|^), one obtains 

1 /--i 
k = —- dK{—X + sinXcosX + XcosXcosKX — sinXcosKX+ 
zvr Ji 

+ KXsinXsinKX) = — (A — smAcosA), (33) 

TT 

which varies from zero to unity as A goes from zero to vr. It is worth noting 
that the middle of the evolution, A = 7r/2, corresponds to a half-integer 
Chern-Simons number, k = 1/2 [0. It is clear now that the boundary 
conditions (^) for the function K in Eq.(|23|) are important. Imposing instead 
the conditions 

K{0) = K{oo) = 1, (34) 

the value of the integral in the right hand side of Eq. (^) would be zero as 
would be the Chern-Simons number. So, the boundary conditions ( plj) insure 
that the field (^) interpolates between homotopically different vacuum sec- 
tors. On the other hand, the conditions (Q) would make the field interpolate 
between vacua inside the same vacuum equivalency class. 

Now, being aware that the field (|23| ) indeed relates to the transition be- 
tween distinct vacua, we adopt a slightly different interpretation of this ex- 
ample. Let us treat the quantity A in Eq.(|25D not as a function of time but 
rather just as a parameter running from zero to vr. This means that the field 
( P5D is viewed now as a one-parameter family of static field configurations or 
a path in the EYM configuration function space. 

Next, we want to specify a gravitational field related to this path. The 
most natural choice is for the matter field ( p3D and the corresponding gravi- 
tational field to form an admissible set of the initial data on an initial time- 
symmetry hypersurface. So we impose the initial value constraints, from 
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which the only non-trivial one in our case is Einstein's equation ([B8|). For 
the sake of convenience only, we impose also the equation (piO|) , which spec- 
ifies the lapse function a. Using Eqs. (|BT|), (^), (|BB|), (|BTI]| ), we get the 
following equations for the spherical gravitational field related to the path 



m 



m' = stn^X{{l )K'^ + stn^X- ^-^}, 

a' = 2sin^X a. (35) 

r 

Straightforward integration of these equation with the boundary conditions 
( p!^ allows us to express m and a in terms of the quantities, which parame- 



terize the field (H) 0: 



air 



m(r] 



exp{-2sin'^X / K'^ — }, 

Jr r 



cr r 



f{K'^ + sin'X ^ ,^ )adr. (36) 



If K{r) is a sufficiently smooth function then these expressions specify a 
metric, which is asymptotically fiat and regular for any A and takes the 
vacuum value, rj^^, at the ends of the path, A = 0, vr. The corresponding 
(dimensionless) ADM mass is 

M = lim m(r). (37) 

r— >oo 

Recall now that K{r) in the above expressions may be arbitrary. This means 
that Eqs.(^),(^) define not one path but rather a family of vacuum-to- 
vacuum paths 

{Aj{x,X);g^,{x,X)}\^^^^, A G [0, vr] (38) 

in the EYM configuration function space, each path of the family being 
specified by a function K{r) satisfying conditions (P^. The energy (^) for 
each path vanishes for the vacuum values, A = 0, tt, and has a maximum in 
between at A = 7r/2. 



Now, one can see that the fields (p3D,(p6D given at the point A = tt/2 



coincide with the solutions of the EYM equations found by Bartnik and 
McKinnon (BK) provided that the function K{r) is identified with the 
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corresponding BK magnetic amplitude w„(r) (see Appendix B). If the num- 
ber n is odd, then K{r) = Wnir) obeys the boundary conditions (|2^), so 
odd-n BK solutions belong to paths connecting distinct vacua (even-n so- 
lutions belong to paths interpolating between representatives of the same 
vacuum class as Wn functions in this case obey the condition (|34D ). This 
circumstance has allowed us initially to interpret the odd-ra BK solutions in 
Ref ■ I T0| as sphalerons, because they relate to the top of the potential bar- 



rier dividing homotopically distinct EYM vacua, A = 7r/2, and also have the 
Chern-Simons index one-half. However, as is mentioned in the Introduction 



(see also |T^), the analogy with the standard YMH sphaleron is incom- 
plete because none of the BK solutions relate to an absolute minimum barrier 
height, and hence an additional investigation is required in order to establish 
the sphaleron nature of these solutions. 



4 The winding number changing transition 
amphtude. 

In this section we represent the amplitude of a quantum transition between 
homotopically non-equivalent EYM vacuum sectors at non-zero energy a la 
Bitar and Chang. We assume that the transition occurs only along the paths 
(|23D , (|36D , which is our basic approximation. Notice that Bitar and Chang's 
approach gives an adequate description of the transition process provided 
one is able to take into account all VTV paths. But the main advantage of 
the method, perhaps, is that it may give a good description of the transition 
even if one takes only some particular set of the vacuum-to-vacuum paths 
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To obtain the transition amplitude, let us first choose a path from the 

family (p3| ) ,(|36| ) . To find the partial amphtude related to this path, allow for 

the parameter A to depend on time: A -^ A(t), and calculate the action of 

the fields. Inserting Eqs.(p3| ) ,(|36| ) with A = A(t) into Eq.(l), represent the 

action in the form 

Atx f 
Seym = ~ / (-^g + Lym) dt, (39) 

where the Lagrangians {47r/g^)LG and [Att / g"^) Ly m of gravitational and the 
YM field are the spatial integrals of the corresponding Lagrangian densities 
given by Eqs.(H),(|^). Consider first the matter part of the action. Using 
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Eqs.(^, (p6D, ( P2| ) and taking into account the Einstein equations (pSj), one 
obtains 

Lym = ^A^ - r m'adr, (40) 



where 

^(A) = / !-{K'' + 2sm2A 7 ^ ^ )^r, (41) 

Jo cr r^ — 2mr 



quantities m, m', a being given by Eqs.(|35D, ( p6D and /('(r) being a function 
specifying the path under consideration. 

The calculation of the gravitational part of the action is performed in the 
Appendix A. For regular geometry described by Eq. ([l3D the horizon term in 
Eq. (|A20|) disappears, so one gets 



^G 







ma'dr (42) 



(when one passes in Eq.( [A20D from dimensional variables t,r,m to dimen- 
sionless ones, which are used throughout in the paper, the multiplier 4:n/g'^ 
explicitly depicted in Eq. (|39|) arises). 



One can see that the integral entering Eq.(J2) combines with that coming 
from Eq.(l^) to give the total derivative under integration. This allows us 
to represent the total EYM action in the following form 

Seym = ^ fi^X' - U{X)) dt, (43) 

g J z 

where the effective mass term, /i(A), is given by Eq.(^T]), and the potential, 
f/(A), reads 

POO 

U{\) = / {ma)'dr = m(oo). (44) 

Jo 

This coincides with the ADM mass given by Eq.([37|) (that is true provided 
that the geometry is everywhere regular; see Appendix A). Explicitly one has 

UiX)^U[Kir),X] = 
= sin'X / {K'^ + sin'X ^ ' )exp{-2sin'\ / K'^—)dr. (45) 

JO ZT Jr T 

One may see that the result obtained corresponds to the action of an effec- 
tive particle with position-dependent mass, /i(A), moving in a one-dimensional 
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external potential U{X). The potential has the typical barrier shape: for each 
K{r) it vanishes for the vacuum values, A = 0, vr, and reaches a maximum in 
between at A = 7r/2 (the latter can be seen if we pass to a new independent 
variable z = r/sinX under the integration in Eq.(^5D). 

Thus, we arrive at the one-dimensional barrier transition problem. The 
corresponding one-dimensional Schrodinger equation then reads 

with p being the momentum conjugated to A, and the quantity E has the 
sense of the energy of the asymptotically free quantum states (the operator 
ordering problem can be avoided in this case |ll[|). This allows us to write 
down the corresponding partial WKB transition amplitude as follows 



4:71 4tt r'^ I 

AKi^r) = Bexp{t—^[K{r)]] = Bexp{%— dA J2/i(A)[E - f/(A)]}, (47) 

Q^ O"^ Jo ^ 



where B absorbs all other WKB factors, which are inessential for our present 
considerations. The subscript K{r) indicates that this partial amplitude 
relates to a single path specified by K{r). If the quantity E—U{X) is negative, 
then one should take the value of the square root lying in the upper half of 
the complex plane. 

The last step in finding the transition amplitude is to represent the total 
amplitude as the sum over all partial amplitudes (|i7|): 

A=J2^Kir). (48) 

K(r) 

This expression gives the amplitude of the winding number changing tran- 
sition in the EYM theory at arbitrary energy E. It is clear, however, that 
taking the sum is virtually impossible. The only way to estimate this sum is 
to make use of the stationary phase approximation. However, as was shown 



in Ref.|TT|, the finding of the exact stationary phase path implies solving a 
system of the coupled differential equations, which, unfortunately, lies be- 
yond our abilities. 

To proceed further we will not calculate the amplitude, but consider in- 
stead a more simple problem. Namely, we want to find out under which 
conditions this amplitude will not be small. 
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We use the following terminology. Let the energy, E, be fixed. Consider 
a path (p3D ,(|36D that lies entirely in the classically allowed region: 



U[K{r),\]<E, AG[0,7r]. (49) 

Such a path will be called an overbarrier path. If a path passes also through 
the classically forbidden region then we shall call it an underbarrier path, for 
the evolution along such a path implies barrier penetration. 

For an underbarrier paths the amplitude (^Tf ) is small, as it includes the 
small tunneling factor. This allows us to exclude from the sum ( ^Hf ) the 
contribution of all underbarrier paths, as this contribution is certainly small. 
The remaining sum over the overbarrier paths will not be small if only this 
sum includes a term (or terms), AK(r), whose value is stationary with respect 
to small variations of K{r). Such a term corresponds to the contribution 
of a stationary phase path. Thus, the amplitude (^) will not be small 
provided there exists a stationary phase path on the set of all overbarrier 
paths. Therefore, to proceed further, we first need to select the overbarrier 
paths from the set (^3]) , (|36|) , and next to look for a stationary phase path 
among them. 

5 Structure of the energy surface 

To select the overbarrier paths in accordance with Eq.(^) one has to know 
properties of the potential U[K{r), \]. It is useful to treat the functional 
U[K{r), A] in geometrical terms, viewing it as an infinite dimensional surface 
in the corresponding function space. Let us call this surface an energy surface 
or a potential barrier surface. It is natural to call the directions in this surface 
generated by d/dX and d/dK{r) a transverse direction and a longitudinal 
direction respectively. Indeed, changing of A with fixed K{r) corresponds to 
motion across the barrier towards a neighbouring vacuum. Changing of K{r) 
implies passing to other paths, i.e. the motion along the barrier. When K{r) 
is fixed, the potential reaches a maximum at A = 7r/2, so we will call the 
functional 

e[K{r)] = U[K{r), A = 7r/2] (50) 

a barrier height functional and say that it defines the profile of the top of the 
barrier. 



Consider the critical points of the energy surface. It is clear that such 
points have to belong to the top, A = 7r/2. Direct variation of the functional 
(BO) yields (see Appendix C) 



Se = 2 /"{-((I - —)aK'y + ^^^, ^K }SKdr, (51) 

Jo r r^ 

where 6K{0) = 6K{oc) = 0, functions m,a are given by Eg. ( |C3| ) . The 
vanishing of this variation requires the condition 

^^^.?I^^,K'Y^a!SlE!f2l (52) 

to hold, which is exactly the YM equation ( |B11| ). Noting then, that the 
functions m, a in this case satisfy the equations 

a' = 2— a, {ma)' = {K" + ^ ^^^ > )cr, (53) 



which are identical to the Einstein equations (P12|) , (|B13| ), we conclude that 



the condition 6e = insures the on-shell condition 6Seym = 0. All points 
of the energy surface relate to the regular field configurations, hence, critical 
points of the top relate to regular solutions of the EYM equations, i.e. to 
the EYM sphalerons. One may say that the EYM sphalerons "lie" on the 
top of the potential barrier. (Of course, these arguments only establish the 
coincidence of the extrema of the action and those of the truncated mass 
functional under consideration. When all field degrees of freedom are involved 
rather than only those which are spherically symmetric, the correspondence 
between extrema of the energy and the action is not proven, although it is 



very likely that it indeed holds; see Ref . [p5|] . ) 

Next, we consider the local properties of the energy surface in the vicinity 
of a critical point. Firstly, we know that the potential energy decreases when 
the transversal A-coordinate deviates from the A = 7r/2 value, so locally, 
near any critical point, the energy surface is a saddle, possessing at least one 
negative direction. Therefore, each EYM sphaleron has at least one negative 
mode which we will refer to as transverse rolling down mode [^. In the YMH 
theory this mode is the single negative mode of the sphaleron (provided that 
the Higgs field self-coupling is not too large ]19[|), so the energy increases 



in any longitudinal direction. In the EYM theory there exist additional 
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sphaleron negative modes, for the potential energy decreases also in some 
longitudinal directions. Indeed, direct calculation (see Appendix C) shows 
that in the vicinity of a sphaleron the following expansion holds 

e[wn{r) + ip{r)] = e[wn{r)] + 6^e + ..., (54) 

where the first order term vanishes and dots denote higher order terms. The 
second order term reads 

^'^ = 1 vi-^ + V)v>dn, (55) 

with the new radial coordinate, r,,, being defined by the following relations 

^ = a(l — 2m/ r), r^{r = 0) = 0, the perturbation '^{r^) obeys (/^(O) = 
^9(00) = 0. The effective potential is 

V = a\l ){^ — + -w'^Wn{wl - 1) - -<2(i _ ^ " ^ )}, (56) 

where functions m and a relate to the corresponding sphaleron solutions, 
the derivatives are calculated with respect to r. One may see that the 
scalar product for different perturbation modes can be naturally defined as 
< V'i,V^2 >= J fi'^2dr^; independent modes being orthogonal. It is obvious 
that if the differential operator in Eq.(^5|) has negative eigenvalues, 

(-^ + V)V = ^V, ^' < 0, (57) 

corresponding eigenfunctions, ip, will specify the longitudinal negative di- 
rections on the energy surface, that is, the longitudinal negative sphaleron 
modes. It is worth noting that Eq.(|57D coincides (up to a constant multi- 
plier) with that first obtained by Straumann and Zhou |2^ in their analysis 
of the linear stability of the BK solutions, provided that one identifies the 
quantity u with the time frequency of perturbations of the background BK 
solution. (Notice, that Straumann and Zhou derived Eq.(0) linearizing the 
EYM equation near a background equilibrium solution. It is interesting, that 
one may obtain the same equation from the energy considerations. A similar 
situation is also encountered in the theory of relativistic stars [Q , where the 



dynamical pulsation equation for a star composed from perfect fiuid can be 
derived by making use of the two different approaches.) 
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It is known ^^, that the equation (p7\ ) has exactly n independent neg- 
ative eigenvalue solutions. Thus in the vicinity of the n-th sphaleron the 
energy surface has, apart from one transversal negative mode, also n longi- 
tudinal negative modes. The total number of the negative modes, n + 1, is 



in agreement with the result by Sudarsky and Wald |25l . 

It is obvious now that when the energy E exceeds an extremal (non-zero) 
value of the barrier height, that is, the mass of the n-th EYM sphaleron, 

E>e[wnl (58) 

then the condition (^91) is satisfied provided that K{r) is close to Wn{r), 
i.e. K{r) = Wn{r) + f{r) with ip being a small perturbation. This means 
that the path passing through the n-th sphaleron (such a path is defined 
by Eqs.(^),(^) with K{r) = w„(r)) as well as neighbouring paths will be 
overbarrier. 



6 The stationary phase path 

Let the condition (^) hold. Consider the overbarrier paths defined in the 
preceding section. Now, we are looking for the stationary phase path on the 
set of these overbarrier paths. Our aim is to show that the path passing 
through the sphaleron will be such a stationary phase path. This means that 
the quantum phase $[i^(r)] defined by Eq.(^7D with K{r) = w„(r) should 
be stationary on variations of K{r). To check the stationarity of the phase 
let us consider an arbitrary variation oi K = Wn'- 

K{r) = Wn{r) + aip{r), (59) 

where a is the variational parameter, and the perturbation obeys ip{r) = 
Olr"^) as r — > 0, ip{r) = 0{l/r) as r ^ oo. Inserting this into $[i^(r)], 
U[K{r-), A] and £:[-ft'(r)] defined by Eqs.(^7D,(^5D and ( |5D| ) respectively, one ob- 
tains the phase $<^(a), (index ip refers to the choice of the perturbation), and 
the one and two-dimensional sections of the energy surface: e^(a;), f/^(a;. A). 
The phase ^[/^(r)] will be stationary if $^(q;) has an extremum for any (f{r) 
in Eq.(|9D. 

One should note that, strictly speaking, the exact stationary phase path 
is not contained in the path family (|23|),(|3^). Finding such a path implies 
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solving the general equations of motion |TT[, which lies beyond the scope 



of our present analysis. Paths (p^),(p6D specified hj K = Wn may be only 
approximately stationary. This means that for these approximate paths and 
any independent perturbations ip{r) in (p9D, the positions of the extremum 
of the corresponding phase functions $^,(0;) do not coincide exactly with the 
position of the sphaleron, a = 0, however they are close together. In other 
words, functional derivatives 6^[K{r)]\j^^^ although do not vanish exactly, 
are nevertheless small. 

Consider first such perturbations ip{r) which increase the barrier height: 
e{a 7^ 0) > b{0); there are infinitely many such perturbations. One may see 
that the phase is stationary with respect to all these perturbations. Indeed, 
in this case the corresponding two-dimensional section of the energy surface, 
U^{a,\), has the typical saddle shape shown in Fig.l. The saddle negative 
A-direction on the picture (shown by the horizontal arrow) specifies the trans- 
verse rolling down mode of the sphaleron |^6[|. It is clear from this picture 
that the path passing through the sphaleron is the minimal potential energy 
path interpolating between distinct vacua, so it corresponds to an extremum 
of the phase. The analogous situation takes place in the YMH theory where 
the sphaleron is the highest-energy point on a minimal-energy path connect- 
ing distinct vacua (at least for small values of the scalar field self-coupling 
[p!9[|), which implies the stationarity of the phase for such a path 0. Direct 
numerical inspection shows that, for several ground state sphaleron energy- 
increasing perturbation modes ip{r) checked, the phase $^(a) indeed has an 
extremum (minimum) in the vicinity of zero value of a. 

However, contrary to the situation in the YMH case, in the EYM theory 
not all perturbations (|59D increase the barrier height, for EYM sphalerons 
possess also longitudinal negative modes. Inserting such perturbation modes 
in Eg. (|59D , one obtains the corresponding two-dimensional sections of the en- 
ergy surface which have the typical shape shown in Fig. 2. One may see from 
this picture that, for such a section, the path passing through the sphaleron 
is not the minimal energy path. So, with respect to these negative pertur- 
bations, the stationarity of the phase is not obvious (it is usually assumed 
that a sphaleron has to have one and only one negative mode in order to 



be significant for the transition processes [|19|). Nevertheless, we are able to 
demonstrate the stationarity of the phase also in this case, at least for the 
ground state (n = 1) sphaleron. 

For the ground state sphaleron there exists only one negative eigenmode 
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to Eg. (|57|) . It turns out that the corresponding perturbation (^) may be 
related to the rescahng of the sphaleron field (the scaling arguments in the 
EYM theory have also been considered in Ref . |5^ . Namely, instead of Eq. ( |5^ ) 
let 

K{r) = wi{pr), (60) 

where (3 is the scaling parameter. Inserting this in Eq.(|^) one obtains the 
following two dimensional section of the energy surface: 

U{(3, A) = ;3sin2 A /"« + sin^ A ^"^^"^^ )exp{-2p^stn^X H w[^—)dr, 

Jo 2iT Jr T 

(61) 
(here Wi = Wi{r)), the value f/(l,7r/2) being the sphaleron mass. This 
function tends to zero when /? — *> 0, oo as well as when A ^ 0, vr, so one 
may see that the potential barrier indeed can be arbitrarily low. The plot of 
this function is depicted in Fig.2. This picture demonstrates explicitly the 
existence of the two independent sphaleron negative modes. One such mode 
is the transverse rolling down mode [^ (the A-direction on the picture), the 
other is the longitudinal negative mode |2^ (the /3-direction). (Note that 
the infinitesimal scaling mode, ip = i9/3Wi(/5r)|„ -^ = rw[{r), is not the exact 
eigenmode for the Eq. (^7D , but rather a superposition of the true negative 
eigenmode and some other mode). 

The scaling behaviour of the quantum phase is defined by the inserting 
into (^): 



^scaUP) = rdX.j2fi{P,\)[E-U{P,\)], (62) 

'J u 

where U[P, A] is given by Eq.(^) and 

M/3,A) = -/ -« + 2.m^A P /^ )dr, (63) 

p Jo a r^ — 2pmr 



with 



air 



P Jo a r^ — 2(3mr 

exp{—2j3'^sin'^\ / w'? — }, 
Jr r 



m(r) = "^ / {w'^ + sin^ y ^ ^ ' )adr. (64) 



Now, one may see that the phase (|62D indeed has an extremum. First note 
that when [3 is small, the phase ^scaieiP) ~ '^1 \fP-, i-e. it diverges when 
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/? — > (remember that the energy E exceeds the sphaleron mass, that is, the 
maximal value of the potential drawn in Fig. 2). Also, and this is the crucial 
point, ^scaieiP) diverges when /3 tends to some value /?c = 1-465, because the 
integral entering Eq.(|63|) diverges in this limit. The value (3 = (3c corresponds 
to such a rescaling of the sphaleron field when the rescaled configuration 
begins to acquire an event horizon (all of the paths (p3D,(|36D,(|60|) with (3 > Pc 
pass through virtual black holes and give an infinite contribution to the 
action). At the horizon the quantity r^ — 2(3mr entering the denominator in 
Eq. (|63|) vanishes, so the integral diverges as does the phase. So, somewhere 
in between, < /? < /?c, the phase must have a minimum (see Fig. 3). 

Thus we can see that when the energy exceeds the ground state EYM 
sphaleron mass, the path passing through the sphaleron insures extremality 
of the quantum phase with respect to any small variations inside the path 
family (PHj) , (pij]) . The existence of this extremum insures that the sum in 
Eq. (|48|) , being evaluated via stationary phase approximation, will be not 
small, as it will include the contribution of a stationary point. Therefore, 
when the energy is large, the winding number changing transition in the 
EYM theory is not suppressed. 

It turns out that the analysis carried out above for the ground state 
sphaleron remains valid for the higher (odd n > 1) sphalerons. This means 
that each higher sphaleron introduces a stationary phase path, only if the cor- 
responding phase will be stationary also with respect to the additional n — 1 
longitudinal higher sphaleron negative modes, which is very plausible. Thus, 
if the energy exceeds the masses of all the EYM sphalerons (i.e. E > 1 m 
the dimensionless units used), the sum in Eq.(H) includes the contributions 
of (infinitely) many stationary points (one point for each higher sphaleron), 
which may lead to an additional enhancement of the transition rate. 

When the energy is less then the ground state sphaleron mass, the situ- 
ation changes. The overbarrier paths in this case also exist, as can be seen 
in Fig. 2. However, there exists no stationary phase path, so the overbarrier 
passages are suppressed by the destructive interference. This is always occurs 
also in the fiat space limit, G — ;> 0, because in this case the threshold energy 
value - the sphaleron mass - being proportional to Planck's mass, goes to 
infinity exceeding any finite value of energy, so that the transition is always 
suppressed. 
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7 Conclusion 

The main results of our analysis may be summarized as follows. In the 
EYM theory, as well as in the pure flat space YM theory, there always exists 
an opportunity to pass over the potential barrier separating homotopically 
distinct vacuum sectors. However, at low energies all the overbarrier histories 
are suppressed by the destructive interference. In the pure YM theory the 
situation remains the same for any energies. In the EYM theory on the other 
hand, when the energy is large and exceeds the ground state sphaleron mass, 
the constructive interference occurs instead. This means that the rate of 
the inclusive ^ fermion number violating reactions at very high energies is 
not small. Such reactions lead to the formation of intermediate sphaleron 
states, the decay of which will produce a large number of gravitons and 
gauge bosons [p^ as well as extra fermions due to the anomaly. It has 



been argued recently by Gibbons and Steif ||3^, that the decay of the BK 
particles should be accompanied by the fermion violation. Our arguments 
show that, at high energies, the probability of BK particles being born and 
subsequently decaying is not small. Notice that the energies available are 
very large - the masses of the EYM sphalerons are of the order of Mpi/g. 
Such processes may arise naturally within the context of superstring theory 
leading to a "primordial" fermion asymmetry. Indeed, lower order terms 
of the expansion of the superstring action in the string tension, give rise to 
the coupled Einstein- Yang-Mills-Dilaton (EYMD) equation of motion. These 
equations possess classical solutions, which resemble in many respects the BK 
solutions and coincide exactly with them in the vanishing dilatonic coupling 
limit |3^ . It is very likely, that these EYMD particles may also be interpreted 



as sphalerons [Q, ^, and our present analysis may be extended to that case. 
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Appendix A The gravitational action 

In this paper we are interested in Lorenzian gravitational action. In view of 
the importance of this issue, we precisely specify in this Appendix our choice 
of action and its relation to the other possible forms of action. 

Our basic requirement for the action is that it is completely free from 
second derivatives of the metric. Otherwise, all the analysis carried in the 



main text becomes cumbersome. As is known |37|, the second derivatives 
entering the Einstein-Hilbert gravitational action, Ry/^, may be combined 
to form a total divergence: 

R./^ = TT./^-d,{V^Wn, (Al) 

where the two-gamma term is 



rr = ^'^^(r2„r" -r-rf,), (A2) 



and 



It is natural then to choose for the gravitational action the following expres- 
sion 

- I6^ /(fl^H- 9„(^W))<(*. = —j^^^S ^^-. (A4) 

because it is manifestly free from all second derivatives. Let the integration 
in this formula be performed over a 4-volume with a boundary S - we are 
working in an asymptotically flat spacetime and assume the S is a distant 
closed boundary which is shifted to spatial infinity in the end of the calcula- 
tions. Then (^1) may be represented in the following equivalent form 



^° = -16^ y ^^''^ - I«3 i 'VW'-VPSM^-. (A5) 

where the integration in the second term is performed over the boundary 
3-surface, whose unit outward normal is denoted by n'^. This expression is 
our choice for the action. The boundary term in (|A5|) is non-covariant, and 



quasi- Cartesian coordinates are implied for calculation of it [0] - we imply 
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that such coordinates exist. To represent this term in a more conventional 
form, assume that there exists an extension of the n^ vector field into an 
open neighbourhood of the boundary; then one may write down 

VX = V^n^ = 9^n^ + T^;X = d^n^ - ^^^r^,n., (A6) 

which allows us to rewrite the surface term in Eg. {\KB^) as follows 

- Y^ fin-^% - g'''K,n,)^\d'x = ^j{K + C)dJ:, (A7) 

where 

K = -V^n^ (A8) 

is the trace of the second fundamental form of the boundary which does not 
depend on the choice of the extension of n^, and 

C=^id^nf^ + d'^n^). (A9) 

is an additional non-covariant quantity. It is easy to see that this latter quan- 
tity does not involve the normal derivatives of the metric. Indeed, represent- 
ing C as C = d^n^ + \g^"rfdag^a^ and decomposing the partial derivative 
operator, (9^, into the normal and tangential components with respect to the 
boundary: 



n„ „„ _„, nun"^ 



d, = ^n^d^ + ((5° - ^V)'9a, (AlO) 



^2 M ^2 



one can see that the normal derivative of the metric in (|A9|) vanishes (the 
condition n^ = const is also to be used). This means that the C-term van- 
ishes on variation of the action ( |X5| ) with the boundary condition Sg^^^lj, = 0, 
so this term does not contribute to the resulting Einstein equations. There- 
fore, if one is interested in the action which just reproduces the correct Ein- 
stein equations on variation, then one may drop the non-covariant C-term 
in Eq.( |A7|) or replace it by the some other quantity giving no contribution 
on variation, for example, by the negative trace of the second fundamental 
form of the boundary imbedded into flat space. In this way one comes to 
the manifestly covariant Gibbons & Hawking's (GH) action |20|. It is clear 



however, that dropping of the surface C-term is equivalent to introducing 
the second (tangential) derivatives into the action, so it is obvious that the 
GH action is not completely free from second derivatives of the metric. 
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Other possible choices of the surface term also do not allow all the sec- 
ond derivatives from the action to be excluded; comprehensive discussion of 
this topic can be found in Ref.|TB|. So, we should inevitably keep the sur- 



face term ( [A^ ) as it stands. The resulting first order action is unique [^ 
Thus, Eq. ([A5|) uniquely defines the gravitational action which is free from 
the second derivatives and gives correct Einstein equations on variation. 

A novel feature arises when an event horizon is present. In this case the 
volume integration in Eq. (|A5|) should be performed only over the external 
with respect to the horizon region, i.e. over the part of the 4-volume en- 
closed by the external boundary, E, and the inner boundary - event horizon 
surface - as well. In obvious notations, the result of such integration may be 
represented as follows: 

iOVrG Jhorizon IoTtG JE 



Sg = -T7-?^I. . {rT,/^-d,{./^Wn)d'x-^^-^ln,W^J\^g\d'x. 



(All) 



Integrating the total divergence, one obtains 



Sg = -T^ f TTV^d'x - -^ / n^W^^\d\ 

iOTTLr J (external region) iOTTLr Jhorizon 

(A12) 
where the surface integral is to be performed over an inner boundary in the 
limit when it tends to the event horizon surface, n^ being the unit outward 
normal to the boundary. One may see that the action ( [A5| ) in this case is not 
equivalent to the manifestly first order expression ([A4|) . Thus, when an event 
horizon is present, the action (|A5D does involve second derivatives. However, 



these derivatives do not influence the variational procedure, as they only give 
the horizon surface contribution to the action, and it is this additional term 
which gives rise to the description of the black hole entropy and other related 
quantities [0, ^ . 

Now we perform an explicit computation of the action ( [A5| ) in the spher- 
ically symmetric case. It is convenient to represent the metric in the form 

ds^ = e'^dt^ - e^dr^ - r\d^^ + sin^§ d^^), (A13) 

where u, A are functions of t, r. Pass to the quasi-Cartesian coordinates, 
x^ = r sin i) cos ip,x'^ = r sin i) sin ip,x^ = r cos ■&. The metric then becomes 

ds"^ = e^df - {5ij + ^^—XiXj)dx'dx^ (A14) 
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I tij i/y j 



). The quantity W^ given by Eg. (|A3|) reads in these coordinates 



W" = (iy°, W') = iXexpi-v), {- - (z/' + -)exp{-X)}-). (A15) 

ty ly ty-' 



Let US specify the integration domain in the volume integral in Eg. (|A^) as 
follows. Let the time, t, vary from —to to to- We consider the case when 
an event horizon is present, whose size, rii, may depend on time such that 
'^i^(^) = T_ff(— ^),''"fl'(±^o) = 0. Such a choice is specified by the situation we 
encounter in the main text, when the time evolution of the metric, whose 
action we want to calculate, starts from the flat space metric at t = —to 
and ends again at the flat space value at t = to passing, in general case, 
through the black hole phase in between. So, the integration domain in 
this case is: t G [— to,to],r e [r// (t) , i?o] , the variables t?, (^ spanning the 
spatial two-sphere, S"^, and the limit i?o ^ oo being implied. The spacetime 
boundary is S = S_|_UStUS_. The spacelike parts of the boundary are S-t = 
{t = ±to,r G [0,i?o];''^ ^ [0,7r],(y9 G [0,27r)}; the corresponding unit normals 



are n(^±)^ = (±exp(z//2), 0); the volume elements are dT,± = Jl'^glcPx = 
exp{\/2)drdQ, where dfl = siri'dd'ddip. The timelike part of the boundary is 
Si = {t G [—to, to], r = i?o, '& & [0, tt], 99 G [0, 27r)}; the corresponding normal 
is n{t)^ = (0, exp(A/2)xYr); the volume element is dT,t = exp{u/2)dtdVL. 

Let us calculate first the surface integral in Eg. (^) . Using Eg. ( |A15|) and 
the definitions introduced, one obtains explicitly 



1 r f°° 

— f W^n^dT, = / drr'^Xexp{- 



t=to 

+ 

t=-to 



(A16) 



r=Ro 



/■*0 Z/ + A, . , 2 ,P — X 

+ / dt{2r exp{ ) — [u r +2r)exp{ 

J~to 2 2 

Direct calculation of the curvature scalar gives the Einstein-Hilbert term 

1 /■ /-to rRo . X — V 

Ry/^d X = dt dr{—dt{Xr exp{ )) + 

"'-to Jrff(t) 2 



AtC J J -to JrH(t) 

+Or\\r V + 2r)exp{ ) — 2r exp{ )) + 

+ r{iy' + X'){exp{'^^)-exp{^^^-))}. (A17) 
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Integrating the total derivatives, adding all these terms together and taking 
the limit Rq ^ oo one obtains 

Sg = - Y^(/ RV^d^x + j( n^WdE) = f_°^ Lcdt, (A18) 

where 

^ " ~4G L "^^^ ^ ){em^—) - exp{^—))dr- 



-— ^{(r z/ + 2r)exp{ — - — ) - 2r exp{ — - — ) + 

4(jr 2 2 

+ rij{t)Xr'^exp{ )} 



(A19) 

r=rH(t) 

Finally, introducing functions m{t,r), cr{t,r) such that exp{v) = (T^(1 — 
2Gm/r), exp{—\) = (1 — 2Gm/r), 2Gm(t,rH(t)) = ruit), one may see 
that the term involving time derivatives in ( |A19D vanishes, and the whole 
expression reduces to a remarkably simple formula 



oo ]^ 

ma'dr a{mr)' 

rH{t) 2 



(A20) 

r=rH{t) 



The two terms in the right hand side of this equation relate to the volume 
and horizon terms in Eq.( |A12D respectively. One may check that for known 
spherically symmetric black hole solutions, the expression ( |A20|) yields a 
value of the gravitational part of the action, which gives rise to regular ther- 
modynamical parameters of a black hole 0], |^ . 



Appendix B Spherically symmetric Einstein- 
Yang-Mills fields. 

In this Appendix we list the full set of the EYM equations in the spherically 
symmetric case. It is convenient to represent the gauge-invariant quantites 
in terms of the variables ilo, fii, /. The non-zero components of the YM field 
energy-momentum tensor then are 



2^2' " ' Aa^^' ' "' r^^' ' '' 2r^ 
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To^ = 2^(//' + ffiofii), (Bl) 

where A = r^ — 2'mr. The component T^ can be obtained from Eq.(|Bl|) by 
changing the sign before the first and the third terms on the right hand side; 
other non-zero components are T^ = T^ = —{Tq + T^). 



The Lagrangian density is 


















g C,YM — — 


itrF^.F'^'^v^: 


= 








^2 


-a.nlif^f 


'^D- 


aA 
r2 


(r+/ 


^nl)-a 


ip- 

2r2 


^^ 


(B2) 


The Pontryagi 


in density reads 




















-tli 


)(1- 


-f) + 


2/(M 


-fa 


o)}. 


(B3) 



Non-zero components of the Chern-Simons current given by Eq.(^ are 

^" = ^^{^l(Pl + P2 - 1) + M + (1 - Pl)P2}, 

K' = -^^{Wo{pI+pI - l)+P2Pi + il-Pi)P2}, (B4) 

The current is not, of course, gauge invariant, but its divergence is, and this 
divergence is just the density ( p3D (up to the ^/—g factor). 
The non-trivial YM equations read 

2 2 

(-(fi'o-fii)y = 2^ffio, (B5) 

a Act 

r^ . Arr 

(-(fi'o-fii))" = 2-/^1^1, (B6) 

r^ A(7 r^ r^ Act 

The independent Einstein equations are 

m' = r'^T^, (B8) 

m = -PT^, (B9) 
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^ = ^(t°-t;). (BIO) 

The equations (P5|)-( pT0D admit non-trivial solutions discovered by Bart- 
nik and McKinnon [^. These solutions correspond to the static purely mag- 
netic case, fio = ^1 = 0- The non-trivial EYM equations then read 

a' = lf%. (B13) 

r 

The regular asymptotically flat solutions to these equations are parameter- 
ized by an integer, n\ the corresponding magnetic function, /, is usually 
denoted by w„. When r runs from zero to infinity, the function w„(r) starts 
from unit value at the origin and, after n oscillations around zero, tends 
asymptotically to (—1)" 0- The metric functions, m{r) and cr(r), mono- 
tonically increase from m(0) = to m{oo) = rrtn and from cr(0) = cr„ < 1 
to o"(oo) = 1 correspondingly. The asymptotic value of the mass function, 
rrin, grows monotonically with increasing n from the value irii = 0.828 to 
rrioo = 1, the physical mass of the solutions is M„ = {\/A7r / g) Mpiuin with 
Mpi being Planck's mass. 



Appendix C Variation of the energy functional 

Consider the barrier height functional defined by Eg . (|50|) : 

s[K{r)] = / (K'^ + ^ . ^ )exp{-2 / K'^-)dr. (CI) 

Jq ZT Jr T 

Let K{r) be a sufficiently smooth function satisfying the following conditions 
K'{r) = 0{r) as r ^ 0; K'{r) = 0{l/r^) as r ^ oo. (C2) 

Define for convenience two new functions 

air) = exp{-2 H K'^ — }, mir) = -^ / V" + ^^ 7 ^^ )cydr, (C3) 
Jr r a(r) Jo r^ 
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they satisfy the following boundary conditions 

cr(oo) = 1, o-(O) 7^ 0; m(oo) < cx), mir) = 0{r^) as r -^ 0, (C4) 
and also, by definition, the following equations 



—a, (ma) = (iT + — 



a' = 2— a, {ma)' = (K" + ^-^^-^)a. (C5) 



Consider small variation 

K{r) -^ K{r) + ip{r), (C6) 

where 

(p{0) = (p{oo) = 0. (C7) 

To preserve the boundary conditions, the variation must also satisfy 

^'(r) = 0{r) as r ^ 0, ip\r) = 0{l/r^) as r ^ oo. (C8) 

Put ( |C6| ) into ( pip and expand the result over ip, one then obtains 

e[K + ip]=e[K]+66 + 6^e + ..., (C9) 

where the first variation is 

6e = 2 r{K'^'a + ^^^['^K y, - 2I{K'^ + i^!-lll!)a}rfr, (CIO) 
JO r^ Ir^ 

and the second variation is 



6h = / {acp^ + a ip^ - 8I{K'ip' + ^^-— ^-ip)a+ 



{K^ - 1)2 



+ (8/2 - 2J){K" + ' ' )cx}dr, (Cll) 



dots in (|C9| ) denote higher order terms and the following new functions have 
been introduced: 

I{r)= TkV^, J{r)= r^"^. (C12) 

Jr r Jr r 
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The boundary conditions (|C2D,([C8D imply that 

/(O) < oo, J(0) < oo, /(oo) = J(oo) = 0. (C13) 

Consider the first variation ( |C1(J| ). Using (|C5|), represent 5e as follows 

5e = 2 {K'if'a + -^—, '-^pa - 2I(may}dr. (C14) 

JO r^ 

Integrating by parts one has 



5e = (2K'aip - 4Ima)C + 2 r{-(K'ayip + a— L L^ + 2m(7l'}dr. 

Jo r^ 

(C15) 

Note, that the boundary terms in this expression vanish. Finding /' from 

(|C12|) and integrating by parts once more, one finally arrives at 

6e = 2 /"{-((I - —)aK'y + ^^^' ~ ^V wr. (C16) 

JO r r"' 

This result agrees with Eq. (^) . One can see that the vanishing of the first 
variation implies the following equation 

((i-!^).^r-. ^'^-^' . (C17) 

Assume now that the first variation vanishes and consider next the second 



variation. Using (|C5|) one obtains 

5^e= / {V^ + or — ^^ -^i[K'^' + ^-— '-^)a+ 

Jo r^ r^ 

+ 8P{may -2J{may}dr. (C18) 

Integrating the fourth term in the integrand by parts, and using Eq.( |C12|) , 
one obtains 

8 / P{maydr = SPma - 16 / Il'madr = 16 / ImaK'ip'—, 
Jo Jo Jo r 

(C19) 
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where the boundary terms vanish. Combining this result with the third term 
in Eg . (10181) , one has 



8 / /{I )aK'ip' + -^—^ '-aip}dr. 

JO r r-^ 



(C20) 



Using ( |C 1 7| ) , ( [CT^ ) , (|C5|) , represent this expression as follows 

2m, ,w X, , „ ^,. 2m 



^ 



+ 



-8 r nil - —)aK'^ydr = - 8/(1 - —)aK' 
Jo r r 

f°° ^., 2m, ^^, , ("^ , 2m, ^^.^ ,dr 

+ 8 / /'(I )aK'^dr = -8 / (1 )aK'V¥' — 

Jo r Jo r r 

2m , , 



dr 



POO 

-4/ (1 
Jo 



)a (fip dr, 



(C21) 



where the boundary terms vanish. The fifth term in 



yields 



poo oo poo poo 

-2 / J{ma)'dr = -2Jma + 2 / J'madr = -2 / maip'2 
Jo Jo Jo 



I dr 



ip 
—2ma—Lp 
r 





°° f^ . if'.. , /■°o,2m ,,, , 

+ 2 / [ma — ) Lpdr = / ( a(f ) ipdr. 

Jo r Jo r 



(C22) 



The first term in ( |C18 



IS 



aLp''^dr 



{cjip'ydr, (C23) 

/o JO 

where the boundary terms are zero. Consider also the following expression 

^""^ ' ^^" (C24) 



0= r^^i-fHy^^ydr. 

Jo r 



Adding the equations ( |C21|) -( |C2^ and introducing also the tortoise coordi- 

(C25) 



nate, t^^, 



dr 2m , 

- — = a(l 

dr^ r ' 



one finally arrives at 



where 



6h = / <^{--r^ + V)ip dr 
Jo dri 



V = a(l-^^){2(Al--))'^'-^ol 



(C26) 



(C27) 



which agrees with Eq. (|56D provided that Eqs.( |C5|) , ( PlTD are taken into ac- 
count. 
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Figure captions 



Fig.l Typical shape of the energy-surface two-dimensional section U^{X, a) 
for energy increasing sphaleron perturbation modes if. The surface forms a 
barrier separating distant EYM vacua (A = 0,7r). Position of the sphaleron 
(a = 0, A = vr/2) is shown by the vertical arrow. 

Fig. 2 Plot of the function U{f3,X) defined by Eg. ([61|) . Vertical arrow 
shows the sphaleron position (/5 = 1, A = vr/2). Horizontal arrows correspond 
to the transverse (A-direction) and the longitudinal (/^-direction) sphaleron 
negative modes. 

Fig. 3 Behaviour of the quantum phase ^scaieiP) defined by Eq.(|62|) with 
E = 1. The path passing through the sphaleron is specified by the value 

(3 = 1. 
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